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ABSTRACT 



The problem of determining the European-style option price in the incomplete market has 
been examined within the framework of stochastic optimization. An analytic method based 
on the discrete dynamic programming equation (Bellman equation) has been developed that 
gives the general formalism for determining the option price and the optimal trading strategy 
(optimal control policy) that reduces total risk inherent in writing the option. The basic 
purpose of the paper is to present an effective algorithm that can be used in practice. 
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1 Introduction 



An essential feature of the currently dominant option pricing theory proposed by Black 
and Scholes is the existence of a dynamic trading strategy in the underlying asset that ex- 
actly replicates the derivative contract payoff [1-4]. However, in general, the market is not 
complete, the contingent claim is not redundant asset and therefore its price cannot be deter- 
mined by the no-arbitrage argument alone. The reasons that give rise to an incompleteness 
of market might be very different, for example, mixed jump-diffusion price process for an 
asset [1,5], stochastic volatility [6], etc. 

In recent years there has been a substantial theoretical efforts to give the pricing formula 
for a derivative security for which an exact replicating portfolio in the underlying asset 
ceases to exist. The typical example involving incompleteness is a model in which the stock 
volatility is a stochastic process. Several approaches to the valuation of the contingent 
claim under random volatility have been suggested in literature [6-11]. Typically the pricing 
formulas involve the unknown and what is more unobservable parameter, so-called market 
price of volatility risk. This fundamental difficulty has led the researches to accept the idea 
of uncertain volatility when all prices for contingent claim are possible within some specified 
range [12-15]. 

An alternative method for the derivative pricing in the incomplete market has been 
proposed in a series of papers by mathematicians Miiller, Folmer, Sondermann, Schweizer 
and Schal [16-20] and by physicists Bouchaud and Sornette [21] (see also [22-25]). The basic 
idea is that the fair price of contingent claim can be found through the risk minimization 
procedure. Different criteria for a measuring the risk inherent in writing an option have 
been suggested including the global and local variance of the cost process [16-18] and the 
variance of the global operator wealth [19,21-23]. We refer to the recent survey paper [20] 
for an exposition of the status of research on the incomplete market involving the stochastic 
volatility and risk-minimization. 

Although significant progress has already been made in the option pricing theory involv- 
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ing the risk minimization procedure, still there exist many open problems including how 
to derive an effective algorithm giving the option price and trading strategy. The purpose 
of this paper is to present such an algorithm that can be used in practice. The aim is to 
show how the problem of option pricing based on the risk-minimization analysis can be re- 
formulated in terms of Maier problem and how the stochastic optimization procedure [26] 
based on the Bellman equation can be implemented to give a reliable numerical technique 
for determining both the derivative price and optimal trading strategy. The application of 
dynamic programming approach to option pricing can be found in [18,28,29]. 

2 Statement of the problem 

We assume a discrete-time iV-period world in which the dynamics of the security price S n is 
governed by the stochastic difference equation 

S n +1 = S n + U = 0, 1, N - 1, (1) 

where £ n is a sequence of random variables whose conditional jump density at time n is 
independent of the asset price S m , m < n, and is given by 

p n ^S n ) = ^P{^<^\ Sn }. (2) 

We assume that at time an investor sells an European-style option with the strike 
price X for Cq and invests this money in a portfolio containing A shares held long partially 
financed by borrowing B in cash. The current value of this portfolio is given by 

V = C = A S - B . (3) 

The investor is interested in constructing the self-financing strategy to hedge the option 
exposure. Since for the incomplete market the exact replication of the option payouts by a 
portfolio of traded securities is not possible, the investor cannot completely neutralize the 
risk inherent in writing the option. Hence the problem is to find such a trading strategy that 
reduces total risk to some intrinsic value[16-25]. 
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To proceed further we need an equation governing the dynamics of the self-financing 
hedged portfolio. First we consider the case of frictionless trading (the effect of transaction 
costs will be also examined in this paper). The value of the portfolio V n at time n may be 
written as 

V n = A n S n -B n , n = 0,l,...,N-l, (4) 

where A n is the number of shares of the underlying asset held long during the time interval 
[n, n+ 1) and B n is the amount of money borrowed. At the beginning of trading period n+1 
just before readjusting the position this portfolio is worth 

K+i = A n S n+1 -(l + r )B n , n = 0, 1, N - 1, (5) 

where r is the interest rate. Therefore the change in the value of the portfolio can be written 

as 

K+i = (1 + r) V n + A n (U -r)S n , n = 0, 1, N - 1. (6) 

Following [19] we propose that the investor's purpose is to maintain a self-financing 
portfolio (4) in a such way that at the expiration date N the terminal value of this portfolio 

V N = (1 + r) Vjv_i + Aat-i (6v-i - r) S N -i (7) 

should be as close as possible to the option payoff 

0x(S N )=max(S N -X,O). (8) 

One way to achieve this purpose is to require that the expectation value of the difference 
between the option value and the value of hedged portfolio at expiration is equal to zero, 
i.e. ~E{0 X (S N ) — V N } = while the variance of this difference E (<SV) — ^/v) 2 } as a 
measure for risk should be minimized by the proper choice of the trading strategy; here E {•} 
denotes expectation with respect to the distributions of £ 0) Ci, ••• ) 6v-i- 
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3 Stochastic optimization 

According the ideas of dynamic programming [26,27], the proper choice of the sequence 
controls A n , B n should involve the information aggregation, i.e. the optimal choice of trading 
strategy at each of N time period should be based on the available information about the 
current values of asset price and hedged portfolio. From a mathematical point of view it 
means that one have to find a sequence of functions (so-called optimal control policy) 

A* n = A* n (S n ,V n ), B* n = B* n (S n ,V n ) n = 0,l,..,JV-l (9) 

that minimize the total risk. In what follows we will use (4) to find an optimal value of B n , 
that is. 

B* (S n , V n ) = A; (S n , V n ) S n - V n , n = 0, 1, N - 1. (10) 

Instead of the problem of minimizing the risk subject to the constraint we consider the 
problem of minimizing the modified risk-function 

R X = E {(6 X (S N ) - V N f + A (9 X (S N ) -V n )}, (11) 

where A is the Lagrange multiplier and the average is made again over all £ n . 

Following the dynamic programming approach, we consider first the last time period and 
proceed backward in time. If at the beginning of the last trading period N — 1 the stock 
price is <SW-i and the value of portfolio is Vn-i , then no matter what happened in the past 
periods, the investor should choose such a trading strategy Ajv-i, B N _i that minimizes the 
risk for the last time period. 

Let us introduce the minimal risk for the last period which is a function of the stock 
price Sn-i and the value of the portfolio Vn-i 

h (Sat-i, V n ^) = min {{6 X (S N ) - V N f + A {0 X (S N ) - V N )\ . (12) 

It follows from (1) and (6) that I can be rewritten as 

lo (Sn-i, Vn-i) = 
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min E ?Ar _ 1 (9 X (<SV_i + £n-iS N -i) - (1 + r) V^r-i - A;v-i (£jv_i - r) SV-i) 2 

+A (9 X (S N -i + frr-iSV-i) - (1 + r) V^-i - Ajv_i (^-1 - r) SV-i) . (13) 

By calculating this function we obtain the optimal value of Ajv_! and thereby the optimal 
trading policy A* N-1 (Sn-i, Vn-i) , B N _ t (Sn-i, Vn-i) for the last period. 

At the beginning of time period N — 2 when the stock price is Sjv-2 and the value of 
portfolio is Vn- 2 the investor should readjust the position in a such way that (Ajy_2, B N - 2 ) 
minimize the risk E^_ 2 {/ {Sn-i, Vn-i)} ■ 

The dynamic programming algorithm takes the form of the recurrence relation 

h (Sn-2, Vjv-2) = 

min E 5jv „ 2 {I (S N -2 + &V-2SV-2, (1 + r) V N - 2 + A N - 2 (6v-2 ~ r) S N - 2 )} . (14) 

Ajv_2 

By calculating I\ (Sn-2, Vn-2) we obtain the optimal function A^„ 2 = A^_ 2 (Sn-2, Vn- 2 ) ■ 
Repeating these arguments we can get the Bellman equation for the period n 

I N - n (S n , V n ) = minE 5n {I N - n -! (S n + £ n S n , (1 +r)V n + A n (f n - r) S n )} . (15) 

An 

The last equation can be rewritten in the form 

In-u (S n , V n ) = min / I N - n -i (S n + £S n , (1 + r) V n + A n (£ — r) S n ) p n (Z,S n )d£. (16) 

The attractive feature of the dynamic programming algorithm is the relative simplicity 
with which the optimal trading policy A* = A; (S n , V n ) ,B* n = B* (S n , V n ) , n = 0, 1, N-1 
can be computed. The basic advantage of general algorithm (15) over functional derivative 
technique [21] is that the original problem (11) is reduced to a sequence of minimization 
problems which of them is much simpler than the original one. 

It might seem that the better choice of control in (15) would be a pair (A n , B n ) giving 
the control policy A* (S n ) , B* (S n ) as the functions of the asset price S n only. However the 
self-financing condition gives rise to the restriction 

A n S n -B n = A n _iS n - (1 + r) B n _, (17) 
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which makes the control problem in terms of the pair (A n , B n ) rather difficult. 

It is clear that the function In (So, Vq) is the minimal risk for the optimal trading strategy 
when the initial value of stock is So and the value of portfolio is V . The initial investment 
required to fund the partially hedged portfolio is nothing else but the price of option Co 
which can be determined from the equation 

dI N (S , C ) 

—dC (18) 

The discussion of when the optimal initial investment Co can be considered as a fair 
option price and problems that might arise from that can be found in [19]. 

4 Transactions costs 

Let us now consider the problem of finding the optimal trading strategy and the option 
price in the presence of transactions costs. We know that the effects of transactions costs 
on the contingent claim pricing might be very complex depending on the size of bid-offer 
spreads, the structure of payoff functions, etc. [4,30-32]. Here we suggest a new algorithm for 
a valuation of option price based on the risk minimization procedure. 

We assume a bid-offer spread in which the investor buys the stock for the offer price 

5 (1 + k) and sells it for the bid price S (1 — k). Again we formulate the problem in terms of 
an investor who sells the European option with payout 6x (Sn) an d who employs the trading 
strategy to hedge the derivative. At time a hedged portfolio is constructed by purchasing 
of A shares at the offer price S (1 + k) and borrowing B in cash at the riskless rate 1 + r 
, so that the ammount of money spent for this portfolio including the effect of transaction 
cost can be written as 

V = A So - Bo + kA S . (19) 

It is assumed here that the investor has no initial position in the underlying asset. The 
investor's purpose is to maintain a dynamic portfolio strategy in a such way that the risk of 
his liability (11) is minimal. 
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Above we have derived a Bellman equation (15) when the asset price and the value of 
portfolio have been chosen as the dynamical variables while the number of shares in portfolio 
has played the role of the control parameter. In the presence of transaction costs it is more 
convenient to make another choice of basic variables and controls. It follows from the self- 
financing condition that the stochastic dynamics of the amount of dollars B n borrowed can 
be written as 

B n+1 = (1 + r) B n + (A n+ i - A n + k\ A n+1 - A n |) S n+1 , n = 0, 1, N - 1. (20) 

Let us introduce a new control parameter Q n such that 

A n+1 = n n , n = 0,l,...,N-l, (21) 

then the dynamical variables B n obeys the stochastic difference equation 

B n+1 = (1 + r) B n + (U a - A n + k | n n - A n |) (1 + Q S n , n = 0, 1, N - 1. (22) 

Now we are in a position to formulate the basic problem. Let us denote by In (S , B , A ) 
the minimal risk that can be achieved by starting from the arbitrary initial state Sq,B , A 

I N (S , B , A ) = n min E{(9 X (S n ) - A N S N + B N f + A (9 X (S N ) - A N S N + B N )\ . 

(22) 

Then the principle of optimality yields the general recurrence relation 

In-u {Sm B n , A n ) = 
minE^ {Jjv_ n _i (S n + £ n S n , (1 + r) B n + (Q n - A n 

(23) 

Clearly, the optimal values of Ao and B$ determining the initial investment Vq as a fair 
option price can be determined by 

dI N (S ,Bo,A ) _ n dI N (S ,B ,A ) _ n 

dB - U ' dA " U - {2A) 
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5 Summary 



To conclude, we have formulated the European-style pricing model for the incomplete market 
in which the risk incurred by selling an option cannot be completely hedged by dynamic 
trading. New effective algorithm based on the discrete dynamic programming equation 
has been presented that gives the option price and optimal trading strategy. The method 
accommodates the effects of transaction costs and can easily be extended to price options 
when volatility is random. The uncertain volatility case can be also treated by the stochastic 
optimization procedure. There are several directions to explore by the method presented 
here. First, one may study the case with imperfect state information regarding the asset 
prices. Also, one can study various adaptive control problems. It should be noted that the 
preliminary work done here might be of big practical importance and therefore merits further 
investigation including he computational aspect of our formalism. 
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